Local and spatially extended sub-Poisson atom number fluctuations in optical lattices 
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We demonstrate that ultracold interacting bosonic atoms in an optical lattice with large on-site 
population show sub-Poissonian on-site and inter-site atom number fluctuations. The experimental 
observations agree with numerical predictions of the truncated Wigner approximation. The corre- 
lations persist in the presence of multi-mode atom dynamics and even over large spatially extended 
samples involving several sites. 
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The preparation of ultracold ensembles with well- 
defined numbers of atoms is advantageous for high- 
precision measurements. In quantum metrology [1] the 
reduced uncertainty in the relative particle number be- 
tween two ensembles is a prerequisite for spin-squeezing 
employed in quantum interfcrometric sensors [2, 3]. In- 
terferometers using ultracold atoms also benefit from a 
well-defined total atom number, e.g., due to constant 
interaction-induced level shifts. 

The basic phenomenon that repulsive interactions can 
lead to reduced atom number fluctuations has been ob- 
served in shallow traps [6] and optical lattices [7-11]. 
Lattice systems have been used for preparation of spin- 
squeezed states for quantum-enhanced interferometry [7] . 
Here we report on the suppression of absolute (on-site) as 
well as relative atom number fluctuations below the limit 
for independent classical particles in a few-site optical 
lattice. We work in the large site-occupation regime that 
is different from the low- filling limit where, e.g., Mott- 
insulator physics can be studied. By site-resolved detec- 
tion we observe reduced fluctuations after adiabatically 
changing the lattice parameters. This method provides a 
much lower limit for the absolute number fluctuations 
than recently developed dissipative approaches [4, 5]. 
Our numerical simulations show that the observed fluctu- 
ations cannot be explained by simple models where only 
one mode per lattice site is considered. A multi-mode 
treatment is necessary and our numerical solutions in- 
dicate the presence of quantum correlations over spatial 
regions of several sites which reveal themselves as spin 
squeezing between the considered regions. 

We load a Bose-Einstein condensate (BEC) of about 
(N) = 5000 87 Rb atoms into a ID optical lattice re- 
sulting in significant population of six sites (not ana- 
lyzing two more sites with a population of less than 
5% of N). Squeezed number fluctuations are observed 
after slowly and linearly ramping up the lattice po- 
tential Vi{x,t) = V(t) sin 2 (irx/d) (the ramp speed is 
dV/dt = 4Hz/ms) from V ~ 430 Hz to V ~ 1500 Hz, 
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FIG. 1. Spatially-resolved atom number fluctuations in an ap- 
proximately symmetric six-site optical lattice after slow ramp 
up of the lattice potential, (a) Single-shot absorption image 
of the atomic cloud with typical mean atom numbers indi- 
cated. The experimentally observed normalized on-site fluc- 
tuations £ 2 (dark) are compared to TWA numerics (light), 
(b) Experimental (dark colors) and numerical (light colors) 
results for the relative atom number fluctuations £jy between 
all site-pairs that increase with pair separation (different col- 
ors). Representative one standard deviation statistical error 
bars (transparent boxes in b) show increased experimental 
uncertainty for the outer wells. For this summarizing figure 
we average all measurements taken for final lattice heights 
larger than 900 Hz. 



with the lattice spacing d ~ 5.7 /im. Some experiments 
were done with final lattice heights below V ~ 430 Hz 
in which case the ramp started from Vb — 250 Hz. Ad- 
ditionally the atoms are held in a cigar-shaped trapping 
potential Vzd{x,U,z) = m[uj 2 x 2 + w\{y 2 + z 2 )]/2 with 
trap frequencies ui± ~ 27r x 427 Hz and u> ~ 27r x 21 Hz. 
By single-site-resolved imaging we detect the local atom 
number rij and its fluctuations An 2 in each site i. 

In Fig. la we plot the observed on-site fluctuations 
normalized to the fluctuations expected for a multino- 
mial distribution. In the central sites we find suppres- 
sion of these fluctuations (£ 2 < 1, defined below) due to 
repulsive interactions. The increase of the variance to- 
wards the outer wells is qualitatively explained by ther- 
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FIG. 2. Evolution of the on-site fluctuations during the lattice 
ramp. Experiment (red), eight-site discrete model at T ~ 
16 nK (blue line), advanced model (TWA) at T ~ 5.5 nK 
(black circles) and T = (inset), (a)-(c) show £ 2 in sites 4-6, 
respectively. Temperatures are chosen for the best agreement 
in all experimental observables shown in this manuscript. 
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mal phonons being dominantly localized in those sites. 
Figure lb shows the spatially extended character of the 
number correlations across the lattice by considering rel- 
ative atom number fluctuations £^ (defined below) be- 
tween all well-pairs. We measure strongly suppressed rel- 
ative fluctuations for all combinations of the central four 
wells. We also show numerical results of non-equilibrium 
dynamics based on the truncated Wigner approximation 
(TWA) [12, 13] (as detailed below) providing good agree- 
ment with the experiment for the central well-pairs, but 
underestimating the outer well contributions. This de- 
viation might be explained by experimental noise and 
the ID nature of the numerical model which cannot ex- 
actly reproduce the physics in the experimentally realized 
crossover regime between ID and 3D. It is important to 
note that simpler models, which neglect the on-site multi- 
mode structure, do not reproduce the experimental data. 

The evolution of the on-site fluctuation with lattice 
height is shown in Fig. 2 for three representative wells. 
Between different experimental runs the integral atom 
number A" varies by A A" ~ 1.6^/ (N). This technical 
noise is compensated by normalizing all data to the mea- 
sured A" in each shot. Due to possible long term drifts of 
the lattice position we calculate pi = (rii/N), the multi- 
nomial probability for well i, for each dataset consisting 
of 25-40 experimental realizations. The on-site fluctua- 
tions Q 2 = (n? ( corr ))/Pi(l ~Pi)N are then obtained from 
the corrected site population ni.( corr ) — rii — piN and 
normalized to the variance for a multinomial distribu- 
tion. The results are averaged over a few datasets such 
that approximately 270 experimental runs contribute to 
each data point and the uncertainty is estimated as the 
standard deviation of the different datasets. Photon 
shot-noise in the detection process adds to the measured 
fluctuations which we subtract from the experimental 
data [7]. With increasing lattice depth the central wells 
show stronger suppression of ( 2 , but the fluctuations in 
the outer wells do not decrease significantly with lattice 
height which can be explained by the localization of the 
lowest phonon modes mainly on the edges. 



FIG. 3. Normalized relative atom number fluctuations 
versus lattice height for experiment (red circles), TWA nu- 
merics at T ~ 5.5 nK (black open circles), the local two-mode 
model at T ~ 10 nK (green curve), eight-site discrete model 
at T ~ 16 nK (blue line), (a)-(c) correspond to the well-pairs 
(34), (45) and (56). At bottom row £jy between larger subsys- 
tems of adjacent sites binned together; (d) sites ({23}{45}); 
(e) sites ({123}{456}); (f) the cumulative contribution of the 
number of continuum HFB phonon modes to the correct value 
of £ir for the pairs (34) (dashed line) and ({123}{456}) (solid 
line) at T ~ 5.5 nK (red) and T = (black). 



In future experiments different lattice sites might be 
used as independent mesoscopic gases with highly sup- 
pressed fluctuations of their total atom number. There- 
fore we analyze the local fluctuations without correct- 
ing for fluctuations of N, the integral atom number, and 
reference the fluctuations to the shot-noise limit of in- 
dividual samples. We measure a four-fold suppression 
of the local fluctuations An 2 /(rii) in the central wells 
as compared to the normalized variance of the integral 
atom number AN 2 /(N). The fluctuations are even sub- 
Poissonian (An?/(nj) ~ — 2.2dB), showing that this ex- 
perimental method can be used to prepare individually 
addressable samples with well-defined atom numbers [4- 
6]. 

To probe spatially extended correlations we analyze 
the relative atom number fluctuations between two lat- 
tice sites and between larger subregions of the lattice 
when the atom numbers of the individual sites are added 
together. In Fig. 3a-c we show the evolution of rela- 
tive atom number fluctuations between adjacent sites, 
normalized to the binomial variance: for the sites 
£n = var(rii - rij)/4pij(l - Pi 3 ){rii + rij), where p i3 = 
(rii/(ni + rij)). The experimental results are limited by 
finite temperature effects as in the case of on-site fluctua- 
tions. Remarkably, are also squeezed between spatial 
regions incorporating more sites (Fig. 3d,e), indicating 
how squeezing can be obtained even over extended lat- 
tice regions. 

To identify the physical processes that are important 
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in the observations we describe the system by different 
models with increasing complexity. The simplest con- 
ceivable model is a local two-mode discrete model H = 
Ecifii - hj) 2 /8 - Ejiajaj + a]a,i)/N + A(h, - hj)/2, as 
used in the analysis of spin squeezing by Esteve et al. [7] . 
Our system is far from the tight-binding regime and the 
parameters are calculated from the ground state wave- 
functions obtained from the 3D Gross-Pitaevskii equa- 
tion (GPE) following the recipe in Ref. [18]. For the 
shallow lattice heights the tunneling timescale is in the 
order of a few tens of ms while its rate of change due to 
the lattice ramp is more than a factor of 10 smaller. Fig- 
ure 3a-c reveals that this simple model can be fitted to the 
experimental observations for the central wells by choos- 
ing the initial temperature T ~ 10 nK. For these wells 
the evolution is adiabatic within this model. However, 
assuming constant T throughout the lattice the fluctua- 
tions in the outer wells are not explained (Fig. 3a-c). 

A more complex model yielding the on-site fluctuations 
is the discrete multi-site model H = — J2<% j> JijW^j + 
ajfii) + ^2j[vjhj + \Ujhj(hj — 1)] for which we calculate 
the Hartree-Fock-Bogoliubov (HFB) modes [17]. The re- 
sults are obtained by including eight sites in which case 
the outermost sites are almost empty - adding more 
sites does not change the results. The parameters Uj,Vj 
are deduced analogously to the local two mode analysis. 
Evaluation of the hopping amplitudes J^- is sensitive to 
numerical accuracy, and thus we take Jy for the central 
sites from the calculated Ej and adjust the remaining 
Jij 's such that the ground state populations match with 
the experiment. Figure 2 shows the on-site fluctuations 
providing the best agreement with the experiment when 
taking the initial temperature as a free parameter and 
assuming adiabatic lattice ramping. The results qual- 
itatively reproduce the level of experimentally observed 
on-site fluctuations across the lattice, but the correspond- 
ing values for £jy are slightly higher than the experiment 
in Fig. 3 and the trend of decreasing fluctuations with 
lattice height is not predicted by this model. 

The failure of the discrete models to reproduce all ex- 
perimental observations indicates the importance of mul- 
timodc effects. To reveal these contributions we describe 
the system by a continuum multi-band ID Hamiltonian 
H = J dx[i'^ (T + V - (i)i> + g 1D ^^^/2], for which we 
first calculate the equilibrium state (T, V, and g\D de- 
note the kinetic energy, potential, and nonlinearity) . We 
find the Bogoliubov theory inaccurate, in which case the 
back-action of the excited-state population on the ground 
state is ignored, and we consequently use the gapless 
HFB formalism where the ground-state and the excited- 
state correlations are solved self-consistently [17]. Specif- 
ically, for Bogoliubov and HFB modes we define the num- 
ber operator in the site i as hi = J i dxifj^(x)tp(x) where 
the integration is over the ith site. The site population 
is then (hi) — No J i dxipo(x) 2 + J.dx (S'(x)6(x)), where 



i>(x) = ipo(x)&o +S(x) with (djdo) = Nq and the excita- 
tions 5(x) = J2j>o [ u j( x )&j — v j ( x )c*j] ■ The atom num- 
ber fluctuations in the ith lattice site Snf = (hf) — (hi) 2 
are evaluated analogously 
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Compared to HFB, the Bogoliubov method overestimates 
both the linear and quadratic terms of Sn 2 [the first and 
the second term in Eq. (1)]. For the shallowest lat- 
tice the difference in 5n 2 is close to 10% (over 100%) 
at at T ~ 5.5 nK (8nK). The quadratic contribution at 
T ~ 5.5 nK is more than two times larger than the lin- 
ear one indicating the importance of the phonon-phonon 
interactions on the number fluctuations even within in- 
dividual sites. The large number of modes required to 
describe the fluctuations emphasizes multi-mode quan- 
tum effects in individual sites as opposed to the single 
mode per site in the tight-binding model. This is visual- 
ized in Fig. 3f where we show the cumulative contribution 
of the number of continuum HFB modes to atom num- 
ber fluctuations, derived from expressions analogous to 
Eq. (1). 

In order to model the non-equilibrium dynamics of the 
experiment in the continuum multi-band theory we em- 
ploy the ID TWA. In the TWA implementation we oth- 
erwise follow Rcf. [13], except that wc develop a pro- 
jection method capable of analyzing the multi-mode ef- 
fects in each site and use the multi-band HFB modes 
in TWA to provide a more accurate description for the 
initial state quantum noise. Although the experiment is 
not strictly ID, the axial modes provide the main con- 
tribution to the observed number fluctuations. Previous 
theoretical studies on atom number fluctuation dynamics 
have either concentrated on double-well systems [14] or 
multi-well systems effectively in the tight-binding limit 
[13, 15, 16], so that significant intra-site multi-mode ef- 
fects were absent. In TWA the quantum field operator ip 
is replaced by an ensemble of stochastic fields ipw satis- 
fying the GPE [13]. The Wigner representation returns 
symmetrically-ordered expectation values that we trans- 
form to normally-ordered ones by constructing an eigen- 
mode basis in each individual site and projecting the field 
operator to the several lowest modes at different lattice 
heights s(t) during the turning-up of the lattice. We de- 
fine the amplitude of the fcth mode of the site i as 

a-i,k{t)= / dx [tpi, k (x,t)]*ip w (x,t) , (2) 

where ipw is the stochastic field and <p^fc is the fcth local- 
ized mode function of the well i. Then the site population 
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reads 

{hi) = 52{&l k &i,k) = E [Kk^-k)w - 1/2] (3) 

k k 

where (• • •) denotes the normally ordered expecta- 
tion value of the quantum operators, and (• • -)w the 
symmetrically-ordered expectation values obtained from 
the TWA simulations. Fluctuations are calculated using 
analogous transformations and the phase in each site is 
defined as a spatial average of the (slowly varying) multi- 
mode field fa ~ arg J dx J2k a i,k^i,k{x). We first let the 
initial state in TWA to evolve in a stationary lattice with 
Vo ~ 430 Hz, in order to damp out excitations due to the 
approximate nature of the initial HFB solution. We then 
turn the lattice up calculating the number squeezing dur- 
ing the ramping and average over remaining small oscil- 
lations in the observables. The black lines in Figs. 2-3 
show the TWA results which reproduce the experimen- 
tally observed behavior for T ~ 5.5(0.5) nK. The TWA 
results are very sensitive to T, however, its absolute value 
can not be expected to match with the temperature in 
the experiment due to the different dimensionality. 

Building on the successful comparison between the- 
ory and experiment we discuss the implications for spin- 
squeezing in this system. A sufficient condition for 
the initial state of quantum-enhanced precision measure- 
ment and many-body entanglement is that the spin- 
squeezing parameter between two subsystems £g = 
N(AS Z ) 2 /({S X ) 2 + {Sy) 2 ) < 1 [19], where the pseudo 
spin operator S z - proportional to the population dif- 
ference in the two subsystems - represents the squeezed 
quadrature and the mean spin is aligned along the xy 
plane. Next to number squeezing which measures 
AS 2 , the coherence ({S x ) 2 + (S y ) 2 ) is an important pa- 
rameter here. For adjacent sites £| ~ £^/(cos0) 2 where 
the coherence is experimentally observable by local inter- 
ference measurements [7] after ballistic expansion [20]. 
The measurements yielded £| ~ — 3.8dB for the cen- 
tral sites. The experimentally measured phase coherence 
is notably lower than (cos0) 2 obtained from the argu- 
ment of the TWA wave function (that incorporates con- 
tributions from the entire momentum distribution), in- 
dicating that spin-squeezing may even be stronger than 
the experimentally inferred value. We numerically sim- 
ulated the effect of the 3D ballistic expansion, but the 
expansion itself cannot explain the lower experimental 
values. At T ~ 5.5 nK and V ~ 1750 Hz we calculate 
£s — — 5.5dB for the site-pairs (3-4) while the value at 
T = is £| ~ — 12.5dB. Numerically we also found no- 
table spin-squeezing between non-nearest-neighbor site- 
pairs and between the extended lattice regions examined 
above. 

In conclusion we prepared sub-shot noise ultracold 
atomic samples with large atom numbers in a few-site op- 
tical lattice. The observations are explained by ID TWA 



simulations. By comparison with simple discrete mod- 
els we find that only the advanced TWA model explains 
all observables consistently revealing the importance of 
multi-mode effects and also providing an additional con- 
firmation of the observations. The temperatures used 
to match the different models to the experiment vary 
significantly, indicating that temperature measurements 
based on fluctuations in this complex system are model 
dependent. The numerical studies show that taking into 
account phonon-phonon interactions is crucial and fur- 
thermore predict that spin-squeezing can be generated 
over extended lattice regions. 
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